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We argue that cosmic strings with high winding numbers generally form in first order gauge 
symmetry breaking phase transitions, and we demonstrate this using computer simulations. These 
strings are heavier than single-winding strings and therefore more easily observable. Their cosmo- 
logical evolution may also be very different. 



I. INTRODUCTION 

Interest in cosmic strings 0, Q has recently been 
resparked by new observational and theoretical findings. 
It now seems that cosmic strings are relatively generic 
prediction of superstring theory which is currently 
thought to be the most promising candidate for a theory 
of quantum gravity. They also appear to be unavoidable 
in grand unification based on field theory, if the same the- 
ory is to explain cosmological inflation as well as unifica- 
tion of elementary particle interactions . Discovery of 
cosmic strings would therefore open up an observational 
window to the very early universe, to extremely high en- 
ergy physics and possibly also to quantum gravity. A 
pair of apparently identical galaxies at redshift z w 0.46 
known as CSL-1 [j| caused a great deal of excitement 
because it seemed to bear the hallmarks of gravitational 
lensing by a string. Even though pictures taken by the 
Hubble Space Telescope showed that it is, in fact, nothing 
but two separate galaxies 0, it led to a large amount of 
theoretical work and significant advances in the field 0. 

In 1980's, cosmic strings were seen as a potential expla- 
nation for the origin of structure in the universe , and 
therefore their observational signatures have been stud- 
ied in detail. The strength of these signatures depends 
on the dimensionless number G/i, in which G is Newton's 
gravitational constant and fi is the tension of the string, 
defined as the energy per unit length. So far, cosmic 
string searches have not found anything, which sets an 
upper bound on the tension. The strongest bounds arise 
from the cosmic microwave background [8j and timing of 
pulsars 0, llOj l ; both give Gn < few x 1(T 7 . The natural 
value of G/i in superstring theory models seems to be sig- 
nificantly lower, around 10~ 10 ...10 -9 0. While these 
strings would be observable with planned gravitational 
wave experiments such as LISA, [llj they are clearly out 
of reach of our current observations. 

However, theories with cosmic strings predict rather 
generally a spectrum of different string tensions. This is 
the case in field theory models in the "Type-I" regime, 
where strings with any integer winding number are sta- 
ble, as well as superstring models in which there are 
stable bound states of D and F-strings known as (p, q)- 
strings 3]. As was mentioned in Ref. (12j . cosmic string 



searches are generally constraining a somehow averaged 
value of the tension. To be specific, let us denote the 
tension of a single- winding string by fj,Q, and the ratio 
of the tension of a given string to this by x = ///Mo- 
Because the amplitude of gravitational waves behaves as 
£l g oc (G^i) 2 / 3 in the relevant range the pulsar bound 

actually constrains x 2 I z G\iq rather than the tension G/io 
of an elementary single-winding string. This means that 
if multiple- winding strings exist, these searches can be 
much more sensitive than they appear. 

The cosmological evolution of a network of strings with 
multiple windings is an open problem, although some 
studies have recently been done. For instance, Ref. 
treats the interactions of strings as local events and ig- 
nores their extended nature. Field theory simulations 
of multiple-winding cosmic string interactions fl3l Xv\ 
have shown that interactions typically lead to lower wind- 
ing numbers. For multiple- winding strings to exist, they 
must have been formed with high winding numbers in 
the first place. The aim of this paper is to show that this 
can indeed happen. 



II. SETUP 

We will focus on cosmic string formation in a finite- 
temperature phase transition in the Abelian Higgs model. 
Realistic grand unified theories should behave qualita- 
tively in the same way, and experience with superstring 
theory models seems to indicate that the same general 
principles would apply also in that case |l5j . The model 
consists of two fields, the complex Higgs scalar (j> and the 
gauge field associated with an U(l) gauge symmetry. 
The Lagrangian is 

C = ™ V" + {D^YD^ - A (0*0 - iff , (1) 

where we have used the field strength tensor — 
dfiA u — d u A u and the covariant derivative D u = d u + 



ieA^. The model has two dimensionless parameters, the 
gauge coupling constant e and the Higgs self-coupling A, 
and the dimensionful scale 77, which corresponds to the 
vacuum expectation value of 0. The corresponding equa- 



2 



tions of motion in Minkowski space are 

-2eLm<f>*D v (j), 



-2A(r</>-?/ 2 )^. 



(2) 



To model cosmic string formation in the early uni- 
verse 0, Il6| . we study our model at a non-zero tem- 
perature. It has two phases: the high-temperature sym- 
metric (Coulomb) phase characterized by a massless pho- 
ton field, and the low-temperature broken (Higgs) phase 
which is essentially a relativistic version of a supercon- 
ducting phase. 

As in superconductors, the "magnetic field" B — WxA 
cannot exist in the broken phase. Instead, it gets con- 
fined into Nielsen-Olesen vortex lines, which are our cos- 
mic strings. The flux $ carried by a string is an integer 
multiple of the flux quantum <5> = 27r/e, and this defines 
the winding number Nw = of the string. At zero 

temperature, strings with Nw > 1 are stable if e 2 /A > 1, 
and at high temperatures one requires e 2 /A > 1.7 [l7j . 
In this so-called Type I regime, the transition is of first 
order. 

We investigate the model in the Type I regime. The 
system is initially in thermal equilibrium in the sym- 
metric phase, and the temperature is then gradually de- 
creased. When the temperature falls below a certain 
value, the symmetric phase becomes metastable. Bub- 
bles of broken phase nucleate, grow and eventually fill 
the whole space. In our analytical calculations, we as- 
sume that the transition is strongly first order so that by 
this time their typical size is much larger than the mi- 
croscopic scales such as the thickness of bubble walls or 
their radius at nucleation. 



III. STRING FORMATION 

It is often assumed that cosmic strings are formed 
by the Kibble mechanism [1|, even though it is strictly 
speaking only valid in the global limit e = 0. According 
to this picture, the complex phase of <f> has a random but 
constant value in each bubble. When two bubbles collide, 
the phase angle takes the shortest way to interpolate be- 
tween the two values. Whenever three bubbles meet, it 
is possible that a hole of symmetric phase is left between 
them. There is then a probability of 1/4 that, when go- 
ing from bubble 1 to bubble 2 to bubble 3, the phase 
angle changes by 2tt. If this happens, a string is formed 
when the hole eventually closes. The number density of 
strings per unit cross-sectional area is therefore roughly 
1/i? 2 , where R is the typical bubble size. Note that only 
strings with Nw — 1 can form. Even Nw — 2 would 
require a simultaneous collision of at least five different 
bubbles at the same point, and since we are assuming 
that the wall thickness is negligible, this never happens. 

Away from the global limit, the Kibble mechanism 
does not fully describe string formation (l6lll8| . Strings 
are also produced by thermal fluctuations of the mag- 
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TABLE I: Parameter values used in the simulations. 

netic field, which become trapped in regions of symmet- 
ric phase liallSj- As we will show, this mechanism can 
form strings with Nw > 1. 

In thermal equilibrium the state of the system is given 
by an ensemble of configurations with a Boltzmannian 
probability distribution, p oc exp(— E/T), where E = 
J d 3 xp is the total energy of the system. In the temporal 
gauge (Aq = 0), the energy density is 



P=~ 2 [E< 



B 2 



D(f> 



+ A(<^-t ? 2 ) 2 , (3) 



where E = — 8qA and B — V x A. Most importantly for 
us, B has approximately Gaussian thermal fluctuations. 

Magnetic flux may get swatted between two colliding 
bubbles |18j . but in this paper we will focus on three- 
bubble collisions, because they can produce higher wind- 
ing numbers . When the three bubbles coalesce leav- 
ing a hole of symmetric phase between them, any flux 
that was there initially will remain trapped in the hole 
because flux lines cannot penetrate the broken phase. 
When the hole closes, the flux trapped in it forms a 
string. 

The winding number of the string is determined by 
the trapped flux, with a possible contribution of ±1 from 
the Kibble mechanism, which we will ignore from now 
on. We must therefore estimate how much flux there was 
in the region before it became trapped. At that time, 
it was still in thermal equilibrium with its surroundings, 
and therefore the flux is determined by the thermal equi- 
librium distribution. We can calculate its typical value 
using the saddle point method. In the symmetric phase, 
where cf> — 0, the energy density is proportional to B 2 . 
This means that the minimal energy i? m i n (<I>) for a con- 
figuration with given magnetic flux <!> through the region 
is proportional to <I> 2 . For dimensional reasons it has to 
be of the form E mul ($) — C$ 2 /i?, where R is the typi- 
cal bubble size and C is a dimensionless number of order 
one. 

The typical value $ rms of the flux, given by the square 
root of the variance of the probability distribution p(<&) oc 
exp(-E min ($)/T), is $ rms w VRT. When the hole be- 
tween the bubbles eventually closes, all this flux turns 
into a string with typical winding number 

N w « $rms/$o « (e/27r)VflT« (e/2 7 r)(u/r) 1/8 T 1 / 2 . 

(4) 

Here we have used dimensional analysis to express R in 
terms of the bubble wall velocity v and the bubble nu- 
cleation rate T as R « (v /T) 1 / 4 . 



FIG. 1: Time evolution of {(f> 2 ) for different cooling rates. The 
time axis has been rescaled by the cooling rate, so that it is 
in one-to-one correspondence with the effective temperature. 



a = 2.5 x 10" 4 




FIG. 2: Bubbles of broken phase at the time when {(j> 2 ) = 
0.33 in the simulations with cooling rates (a) a — 2.5 x 10~ 4 
and (b) a = 2.5 x 10 -5 . With lower cooling rate, the bubbles 
are larger and there are fewer of them, because the nucleation 
rate was lower at the time of their nucleation. 



Because R is much greater than the thermal wave- 
length 1/T, Nw can be large. If we, for instance, assume 
that R is comparable to the Hubble radius 1/H, we find 



N w 



2ng 



1/4 



T 



27r( ff *G A/ o) 1 /4' 



(5) 



where g* is the effective number of thermal degrees of 
freedom. If the tension [1q of a single-winding string is 
low, the typical winding number can therefore be high. In 
particular, this means that the typical observed tension 
/i f» NwHo scales as (G/io) 3 / 4 rather than G/io, making 
the strings more easily detectable. For e = 1 and <?* « 
100, this would strengthen the pulsar bound by about a 
factor of two. 



IV. SIMULATIONS 

To test this theoretical picture, we carried out numer- 
ical simulations of first-order phase transitions in the 
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FIG. 3: String networks formed in simulations with 
two different cooling rates a. In both cases, we have 
used a short period of gradient flow to remove short- 
distance thermal noise, and plotted isosurfaces of \(j>\. Wc 
have made movies of the full time evolution available at 
http://www.imperial.ac.Uk/people/a.rajantie/research/heavystrings/ 



Abelian Higgs model. We discretized the theory on a 
256 3 lattice in the standard way [2(j usin g, th e non- 
compact formulation (see, for instance, Ref. 21] ). The 
values of all parameters are given in Table|U We first gen- 
erated thermal initial conditions by evolving the system 
in time using the classical equations of motion and 
generating random values at random intervals for those 
components of the time derivatives 4> and E = — 8qA that 
do not appear in Gauss's law V • E = 2elm.<fi*(f>. 

After 512 units of time and 64 randomization steps, we 
were satisfied that the state of the system was in equilib- 
rium. After that, we added a small damping term a to 
the equations of motion @, and evolved the system in 
time without any random noise. The damping term was 
implemented in such a way that it only affected the com- 
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ponents of the field that do not appear in Gauss's law. ft 
cools the system down as T e g oc exp(— at) and serves two 
purposes: It causes the phase transition to start earlier 
because the bubble nucleation rate depends on the tem- 
perature, and it removes the latent heat released by the 
growing bubbles. Without damping, the system would 
end up in a state of two coexisting phases. We ran the 
simulation with two different cooling rates a — 2.5 x 1CP 4 
and 2.5 x 1CT 5 . 

The time evolution of the variance (</> 2 ) of the scalar 
field is shown in Fig.^ an d consists of three stages: First, 
the system is in the symmetric phase, with a small and 
almost constant (<fr 2 ). Then, at around t « 0.05/cr, there 
is a rapid transition to the broken phase, after which 
the system equilibrates, apart from a small number of 
strings that were formed in the transition. When the 
cooling is slower, the transition takes places at a higher 
temperature. Because the nucleation rate is then lower, 
fewer bubbles are nucleated and they can become larger. 
To confirm this, we have plotted in Fig.[3the regions that 
are in the broken phase at the instant when (cf> 2 } = 0.33. 

Pictures of the string networks produced in the two 
runs are shown in Fig. [3] As expected, slower cooling pro- 
duced fewer strings. The ratio of the total string length 
between the runs is 2.3, indicating a ratio of typical bub- 
ble sizes of 1.5. 

In both cases, one can see several junctions where two 
strings merge into an heavier winding-2 string. This 
would not happen in the Kibble mechanism. The length 
ratio of windings 2 and 1 is higher in the slower transition 
(0.18 vs 0.15), supporting our theory. However, because 
of the low statistics and finite-size effects, this cannot be 
directly compared with the prediction in Eq. which 
furthermore assumes that Nyy 3> 1- This condition is 
not satisfied in this simulation or indeed in any similar 



simulation unless one can use several orders of magnitude 
larger lattice sizes. 



V. CONCLUSIONS 

We have shown that thermal first-order phase tran- 
sitions can produce strings with high winding numbers 
and predicted the typical winding in Eq. Our nu- 
merical simulations confirm the formation of strings with 
Nw > 1. Quantitative tests are not possible with sim- 
ulations of this type because of the computing power 
they would require, but other numerical tests will be re- 
ported in a future publication |l8j . It is not known how 
a string network with different winding numbers evolves, 
but if the heavy strings survive until today, they will 
be much more easily observable than ordinary single- 
winding strings. The current observational upper bounds 
would therefore become much stronger. This will, how- 
ever, require a large scale simulation of the later evolution 
of the string network and is beyond the scope of this pa- 
per. It would also be interesting to extend this work to 
brane collisions in superstring theory models to find out 
whether heavy strings are also produced in that case by 
a similar mechanism. 
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